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Abstract 



The jamming scenario of disordered media, formulated about 10 years ago, has in 
recent years been advanced by analyzing model systems of granular media. This has 
led to various new concepts that are increasingly being explored in in a variety of 
systems. This chapter contains an introductory review of these recent developments 
and provides an outlook on their applicability to different physical systems and on 
future directions. The first part of the paper is devoted to an overview of the findings for 
model systems of frictionless spheres, focussing on the excess of low-frequency modes 
as the jamming point is approached. Particular attention is paid to a discussion of the 
cross-over frequency and length scales that govern this approach. We then discuss the 
effects of particle asphericity and static friction, the applicability to bubble models for 
wet foams in which the friction is dynamic, the dynamical arrest in colloids, and the 
implications for molecular glasses. 

0.1 Introduction 

Just over ten years ago it was proposed ( |Liu and Nagel, 1998[ ) to approach disordered 
condensed-matter systems in a more unified way than is usually done, by starting from 
the observation that many systems — not just molecular glasses, but also many soft- 
matter systems like granular media, colloids, pastes, emulsions, and foams — exhibit 
a stiff solid phase at high density, provided the temperature and the external forces or 
stresses are small enough. This proposal led to the idea of a jamming phase diagram 
for particulate systems, redrawn in Fig. 10.11 

The very idea that there would be a generic jamming phase diagram is less trivial 
than may appear at first sight, since granular media, emulsions and foams are in fact 
athermal systems — they are represented by points in the ground plane of Fig. 10.11 as 
their relevant interaction energies are orders of magnitude larger than the thermal en- 
ergy scale kgT. While thermodynamics generally guarantees the existence of a unique 
phase for molecular systems at equilibrium, athermal systems strictly speaking lack 
the averaging needed to be able to define a unique phase: their disordered jammed 
states are often history dependent, a phenomenon they share with glasses. 

Nevertheless, it has become clear that it is extremely useful to approach these 
systems with a unifying common focus in mind. For particulate systems like granular 
media, the constituent particles essentially have strong repulsive interactions of finite 
range: particles that are in contact do interact and repel each other, those that are 
not in contact, don't. A detailed analysis of the simplest model of this type, that of 
N repulsive frictionless spheres (or discs in two dimensions) in the large- iV limit, has 
revealed that the point at which they just jam into a rigid phase does indeed have 
properties reminiscent of an asymptotic or critical point that organizes the behavior in 
its neighborhood. This is the jamming point, marked J in Fig. 10.11 It is the aim of this 
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Fig. 0.1 (a) The jamming phase diagram; inside the shaded region, for high density, low 
temperature and small shear stress, particulate systems are jammed, i.e., form a disordered 
solid phase with a finite resistance to shear. The point J is the so-called "jamming point" for 
spherical particles with finite-range repulsive forces. It acts like a critical point and organizes 
the behavior in its neighborhood, (b) The type of forces considered in models of frictionless 
spheres: the force f{r) is of finite range, i.e. vanishes when the distance r between the centers 
of the spheres exceeds some well-defined value. The lower case, in which the force increases 
linearly with the compression S, is often referred to as one-sided harmonic springs. The 
Hertzian force of 3c? elastic spheres increases as 5^'''^ and has a behavior as in the upper 
figure. 

chapter to provide an introductory review of these recent findings, and to provide an 
outlook on the relevance of our present understanding for other systems. We aim our 
review of jamming on those aspects that are of most interest from the point of view 
of the main theme of this book, glasses. Hence we focus in particular on the various 
length scales that emerge near the jamming point, on the vibrational properties of 
the marginally jammed state above the jamming transition, and on the connection of 
these findings with the dynamics. 

This approach bears close similarities with investigations of the rigidity of network 
glasses developed in the 1980's. Instead of density and shear stress, Phillips argued that 
the key control variable is the glass composition ( [Phillips, 1979[ ) . Based on a model of 
rigidity percolation proposed by Thorpe, where springs are randomly deposited on a 
lattice ( |He and Thorpe, 1985[ [Boolchand et ai, 2005 ), it was argued that glass prop- 
erties are controlled by a critical point at zero temperature. The jamming transition 
can be viewed as a special case of the generic rigidity percolation problem in which the 
system self-organized to avoid large fluctuations in its structure. Thus the jamming of 
particles is qualitatively different from springs randomly placed on a lattice. In some 
respects, the jamming transition is simpler, allowing the recent conceptual progress 
reported here. 

Unfortunately because of the focus of this book and space limitations, this chap- 
ter cannot do justice to the wonderful recent experimental developments: many new 
systems have become available or have been analyzed in new ways, and are increas- 
ingly driving new directions in the field of jamming. We refer for an overview of recent 
experiments on colloids to the chapter by Cipelletti and Weeks and for a review of 
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experiments on grains and foams to the chapter by Dauchot, Durian and van Hecke. 
We will point to relevant experiments where appropriate, without going into details. 

0.2 Overview of recent result on jamming of frictionless sphere 
packings 

Close to the jamming point, packings of compressible frictionless spheres exhibit var- 
ious anomalies including a strongly enhanced density of states at low frequencies. Of 
course, such packings are a very idealized model system for granular media. We post- 
pone to section 10.31 a discussion of generalizations of this model such as the role of 
particle anisotropy and friction. 

In computer models of compressible spheres, one usually uses force laws of the type 
sketched in Fig. lO.lf b): two spheres of radii Ri and Rj experience zero force if the 
distance Vij between their centers is larger than Ri + Rj , and have a force that rapidly 
increases with the overlap 6ij = Ri + Rj — Tij if < Ri + Rj. A particularly convenient 
and common choice for computer models is the one-sided harmonic spring model, for 
which the repulsive force fij for particles in contact increases linearly with the overlap 
5ij. Another choice that one often encounters in the literature for compression of 

3 /2 

elastic balls, is the Hertzian force law fij ~ S^j . In studies of three-dimensional 
random packings, it usually suffices to take spheres whose radii are all the same, as 
one easily ends up with a random packing, but in two-dimensions one needs to take a 
polydisperse or bi-disperse distribution of discs in order to avoid crystallization. There 
are various ways to prepare static packings of spheres in mechanical equilibrium. One 
method is to place particles randomly in a box at a fixed density and quench the 
system to its closest potential-energy minimum via conjugate-gradient or steepest 
descent algorithms ( [O'Hern et ai, 2003D . Other methods employ a slow inflation of all 
the radio to first generate a packing at the jamming threshold where the majority 
of the particles experience minute forces that balance on each particle. This jamming 
threshold marks the onset of a nonzero pressure and potential energy. After this initial 
preparation the packing is compressed (or dilated) while continuously minimizing the 
energy dO'Hern et ai, 2003] ) ■ Other methods slowly adjust the radii so as to steer the 
pressure in the packing to a prescribed value ( [Shundyak et ai, 2007[ ). We refer to the 
literature for details. 

We note that the onset packing fraction of the jamming transition for an en- 
semble of states depends on the ensemble. For systems equilibrated at infinite tem- 
perature and quenched to T — 0, the onset packing fraction corresponds to that 
of an ensemble in which each state, or local energy minimum, is weighted by the 
volume of its energy basin. In that case, 0c ~ 0.64 for monodisperse spheres in 



three dimensions (O'Hern et al, 2003). However, if the ensemble is prepared by a 
quench from a system equilibrated at a low temperature, the onset packing fraction is 
higher ( [Berthier and Witten, 2009a| [Chaudhuri et ai, 2010 ) . It is clear from the dis- 



tribution of jamming onsets that this must be the case ( O'Hern et ai, 2003[ ). In the 



infinitc-system-sizc limit, the distribution approaches a delta-function at </) w 0.64 but 

^This is somewhat like the Lubachevsky-Stillinger method to generate random close packings of 
hard spheres l |Lubachevsky and Stillinger, 1990[ |. 
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Fig. 0.2 The density of states of vibrational modes of tliree-dimensional soft-spliere paclc- 
ings witli one-sided iiarmonic forces at varous densities: (a) significantly compressed sam- 
ples; (b) close to the jamming density cfjc ( |0'Hern et al., 2003[ ); (c) on logarithmic scale 
( |Silbert et al, 2005[ ). 



with non- vanishing tails both on the high-density and low-density side ( O'Hern et al., 2003 1 
In systems equihbrated at lower temperatures, lower energy states are weighted more 
heavily. Such states have higher values of t^c, since the potential energy increases with 
(j) — 4>c- Thus, the average value of 0c for the ensemble must increase as the system 
is quenched from lower temperatures, in accordance with the results of Chaudhuri, et 



al. (Chaudhuri et al., 2010), and the system can be described as having a line of jam- 
ming transitions extending upwards from </> « 0.64, as in Ref. (Mari et al., 20091. In 
all cases, however, the onset pressure is zero, so the onset is sharp in terms of pressure 
but not packing fraction, and the vibrational properties, etc. are identical to those 
described below. 



0.2.1 The vibrational density of states of packings 

An important concept in condensed- matter systems is the density of states, D{lj). 
D(uj)Auj is proportional to the number of states with frequency between uj and oj + Aui. 
Here uj refers to the frequency of the vibrational normal modes of the constituent 
particles. The density of states concept is also often used for for electronic states, but 
here we focus on the vibrational states of condensed matter systems. For a crystal, 
or in fact for any elastic medium, D{uj) increases at low-frequencies as w''"^, where 
d is the dimension of space. This generic behavior arises from the fact that sound 
modes have a dispersion relation ijj(k) which is linear in the wavenumber k, together 
with phase-space arguments for the number of modes with wavenumber between k 
and k + dfc. The Debye scaling law D{u}) ~ for d = 3 underlies the ubiquitous 
low-temperature specific heat of three-dimensional solids. It serves as an important 
reference for identifying anomalous behavior — e.g., the well-known enhancement of 
the specific heat over the Debye law is an indication in glasses for an excess density of 
excitations at low frequencies. 

To analyze the vibrational modes, one obtains the so-called dynamical matrix famil- 
iar from solid state physics. This dynamical matrix is essentially the second derivative 
of the inter-particle potential, and hence has only non-zero elements for particles that 
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are in contact in the packing (for the one-sided harmonic forces these terms are espe- 
cially simple as the potential is quadratic in the separation). From the diagonalization 
of the dynamical matrix one then obtains the eigenmodes and their eigenfrequencics 
uj, and hence Diuj). 

Fig. 10.21 shows one of the early results for a three-dimensional packing with one- 
sided harmonic forces dO'Hern et ai, 2003{ [Silbert et ai, 2005] ) . Panel (a) shows that 
sufficiently far above jamming, -D(w) vanishes at small frequencies, in qualitative agree- 
ment with the Debye scenario, but as the density (j) is decreased towards the jamming 
density (f>c, the weight at small frequencies increases. Indeed, closer to the jamming 
density D{uj) develops a plateau at small frequencies, see panel (b). When plotted 
on a log-log scale, as in panel (c), the crossover and the emergence of a plateau is 
even more clear: in these data one observes the Debye scaling only for the largest 
densities, while as the packings are decompressed towards the jamming density, the 
plateau extends to lower and lower frequencies. Clearly, the closer the packings are 
to the jamming threshold, the more D(uj) is enhanced at low frequencies and the 
larger are the deviations from the usual Debye behavior. We will discuss in section 
I0.2.3l the cross-over frequency oj* , that separates the plateau from the lower-frequency 
downturn. 

0.2.2 Isostaticity and marginally connected solids 

How do the excess low-frequency modes arise? The answer is related to the fact that 
at the jamming threshold, packings of frictionless discs and spheres are isostatic, 



i.e., they are marginal solids that can just maintain their stability ( 


Alexander, 1998 


Moukarzel, 1998| Head et a/., 2001||Tkachenko and Witten, 2000||O' 


Hern et al, 2003 


Roux, 20001 Wyart, 20051. The origin of this is the following. Consider a disordered 



contact number. Force balance on each particle implies that the vector sum of the 
forces adds up to zero; therefore the requirement of force balance on all particles with 
contacts gives dNc conditions. If we view the contact forces between the particles as 
the degrees of freedom that we have available to satisfy these requirements, then there 
are ZNc/2 such force degrees of freedom. Clearly, assuming no special degeneracies, 
we arrive at the condition 



for a stable packing: Z > 2d. (0-1) 

in order that force balance can be maintained. 

There is a second constraint in the limit that one approaches the jamming point. 
As all forces are purely repulsive, in this limit almost all the individual contact forces 
must approach zero. For force laws like those sketched in Fig. lO.lf b). this mean that 
as one approaches the jamming point, all nontrivial ZNc/2 contacts must obey the 
"just- touching" conditions Vij = Ri + Rj. The number of degrees of freedom associated 
with the positions of the centers of the particles with contact is dNc, so if we think 

■^There is a subtlety here: in a typical packing there is generally a small fraction of "rattlers" 
or "floaters", particles in a large enough cage of other particles that, in the absence of gravity, can 
float freely without any contact. These should be left out from the counting below, and from the 
determination of the average contact numer. 
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Fig. 0.3 (a) Z scales as |A(^|°-^, from ( |0'Hern et al, 2003] ). These results extend those of an 
earlier numerical study ( [Durian, 1995p . (b) ui* scales as lAt/il"'^. From ( [Silbert et al, 2005^ . 
Together these two sets of data are consistent with uj* ~ AZ. 

of putting the particles with contacts in the right place to allow them to obey the 
just-touching conditions, we need to have 



just-touching conditions at jamming: Z < 2d. 



(0.2) 



Clearly, the value Ziso = 2d, the isostatic value, has a special significance. Indeed, the 
above two conditions imply that as packings approach the jamming point, point J of 
Fig. lO.lf a). from the jammed side, e.g., by decompressing them, one will have 



upon approaching point J: Z I. Ziso = 2(i. 



(0.3) 



Note that these results arc independent of the presence of polydispcrsity and the details 
of the repulsive force law, provided it is of finite range and continuous H 

The above argument will be further justified in section lD.2.4l for frictionless spheres. 
Fig. iO.Sf a) shows numerical simulation results for AZ ~ Z — Ziso, plotted on a log- 
log scale as a function of the distance from jamming, A(j) = — 4>c, in both 2 and 3 
dimensions. In both cases AZ goes to zero at jamming: AZ ^ \/ A<j> ( [Durian, 1995[ 
jO'Hern et al., 2002[ jO'Hern et al., 2003| ) . A recent experiment aimed at testing this 
scaling, albeit in a system with friction, can be found in ( [Majmudar et al., 2007| ). The 
isostaticity concept will be re-examinated in section [0.3. II for ellipses. 

Note that according to (|0.1|) the packings at isostaticity have just enough contacts 
to maintain stability. In this sense, they are marginal packings, packings at the edge of 
stability. Moreover, if one imagines a packing with M fewer contacts than dictated by 
the isostaticity condition, one can according to (|0.2[) deform the packing in M different 
directions in the space of coordinates, while respecting the just-touching constraints (so 

^In passi ng, we note that 4>c, the density of a packing of monodispersc spheres (all the same radii), 
approaches l |0'Hern et al., 2003[ l the random close packing density of hard spheres. There is an active 
line of research aimed at analyzing and relating the concepts of jamming and the maximally random 
jammed state, which has been proposed to replace the concept of the random close packed state 
| |Torquato et al, 2000| [Torquato and StiUinger, 2001| |. 
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Fig. 0.4 Three stages of creating a low-energy Ansatz for a vibrational eigenmode of the 
system, used in the argument to derive the flatness of the D{lj) close to point J, as inspired by 
( |Wyart et al, 2005Q[|Wyart et al, 2005 &[ ). In (c), the solid line illustrates the behavior of the 
smooth amplitude with which the underlying floppy mode is distorted, and which vanishes 
at the boundaries of the box. See text for details. 

that particles arc not pressed into each other) . These therefore correspond to M zero- 
energy deformation modes in which particles slide past one another. These modes, 



which are in general global modes, have been called floppy modes (Thorpe, 1983 



Alexander, 1998) 



0.2.3 The plateau in the density of states and the cross-over frequency, 
Lu*, and length scale, £*. 

The development of a plateau in the density of states is intimately connected with the 
approach to the isostatic jamming point. To see this, we argue as follows ( |Wyart et al., 2005a 
|Wyart et al., 20056| . Let us start from an isostatic packing at jamming, sketched in 
Fig. 10.4( a). Now imagine we disregard ("cut") for a moment the bonds across a square 
or cube of linear size £, as sketched in panel (b). There are of order £'^~^ of these 
surface bonds, and since the orginal packing was isostatic, by cutting the bonds at 
the surface we create of order f'^"^ floppy zero-energy modes within the cube. Each of 
these modes will be very complicated and disordered, but that does not matter here. 
Next, as indicated in panel (c), we use each of these floppy modes to create a varia- 
tional Ansatz for a vibrational eigenmode as follows: we take each of the zero-energy 
floppy modes created by cutting bonds, and distort it by a smoothly-varying sine-like 
amplitude which vanishes at the boundaries of the box, right at the point where we 
have artifically cut the bonds. This long-wavelength distorted mode is a good Ansatz 
for a low-energy eigenmode: when we analyze the potential energy associated with this 
Ansatz-mode, the underlying floppy mode does not contribute to the energy — only 
the fact that the mode has been made imperfect by the elastic distortion on the scale 
£ contributes at every bond a distortion energy of order l/£^ (nearby displacements 
differ from that in the underlying floppy modes by an amount of order the gradient of 
the amplitude, so the energy, which involves terms of order the average compression 
squared, is of order l/£'^). Hence this Ansatz mode will have a frequency uj£ = O {!/£). 
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Of course, in a variational calculation, each of these Ansatz modes will acquire 
a lower energy upon relaxation, but the number of them will not change (rigorously 
speaking this argument yields a lower bound on the density of states). In what follows 
we shall call the modes obtained from the distortion of floppy modes 'anomalous', 
as their nature is very different from that of plane waves. Assuming that the energy 
does not shift dramatically, we have created ~ £'^~^ modes with frequency up to 
LU£ ^ in a box with of order ^ i'^ particles. Thus we have 

f,.Di.)^^^]. (0.4) 

If we assume that D{uj) scales as w*? for small frequencies, we get 

which immediately yields a = 0: D{uj) is flat at low frequencies. 

When the packings are compressed so that they have an excess number of bonds 
AZ = Z — Zc, the above line of reasoning can also be followed to obtain the cross- 
over frequency cj* and length scale £* that separate the range dominated by isostaticity 
effects from the usual elastic behavior. Instead of starting from an isostatic packing, we 
retrace the above construction for a packing with given excess contact number AZ per 
particle. The total number of excess bonds within the box of Fig. 10.4( b) then scales like 
AZ i''^. If the number of degrees of freedom Ne ~ ^''"^ created by cutting bonds at the 
boundary of the box is less than the excess number of bonds in the bulk, no zero-energy 
modes will be created. This fact was already noticed by Tkachenko and collaborators 
in inspiring works ( |Head et ai, 200 1| [Tkachenko and Witten, 2000| . Hence we expect 
a cross-over scale £* when the two terms balance, i.e., when 

Ni, ~ {tf-^^ ~ AZ {tf =^ r - l/AZ. (0.6) 

Likewise, according to the argument above, one expects the crossover frequency w* in 
the density of states above jamming will scale as 

a;*~l/r-AZ. (0.7) 

The scaling of the cross-over frequency, uj* ~ AZ, has been well documented from 
a detailed analysis of the density of states, as shown in Fig. 10.21 At finite compression, 
u* is the frequency above which there is a plateau in D[u), and below which Diui) 
decreases with decreasing frequency. Early data obtained as a function of excess density 
A(j) are shown in Fig. 10.3( b). and are found to scale as lo* ~ |A0|°-^. Since AZ ^ 
|A(/)|°-^, see Fig. 10.3( a). this is consistent with the scaling lo* ^ AZ. An explicit plot 
of w* versus AZ can be found in ( [Wyart et al, 2005a| . 

We stress that in the above argument, it is implicity assumed that the bond strength 
k remains unchanged. This is true for one-sided harmonic forces, but not for Hertzian 
interactions which weaken upon approaching the jamming point: / '-^ 5'^^'^ so that 
k ^ 5^^"^ ^ \J A(\). Since frequencies scale as ^fk, the above arguments generalize to 
non- harmonic forces if formulated in terms of scaled frequencies Cj = uj/Vk. 



Overview of recent result on jamming of frictionless sphere packings xiii 




ip=10"'^ 









































Fig. 0.5 Illustration of the fact that the length scale £* increases as the jamming threshold 
is approached. This figure illustrates the response to a force loading in the center, with 
blue (red) lines indicating an increase (decrease) of the force at a contact; the thickness of 
the lines is proportional to the size of the change in force. The left panel is for a pressure 
p = 10~^, the right one is close to the jamming point, i.e., at pressure p — 10"^. In the 
latter case, the fluctuations are larger and extend over a larger region. A detailed analysis 
shows that this range grows proportionally to the length scale £* . From EUenbroek et al. 
pUenbroek et ai, 2006) |Ellenbroek et al, 2009at . 

While the cross-over frequency is relatively easy to extract from the vibrational 
density of states, the cross-over length £* is more difficult to extract from a mode or 
from response data. Physically, i* is the length scale on which the lowest-frequency 
anomalous modes probe the microscopic structure of the solid, and are therefore sen- 
sitive to its fluctuations. Thus one may expect to observe £* in the fluctuations of the 
linear response of the solid, rather than in its mean behavior. Another difficulty lies 
in the fact that the floppy modes which form the basis for the behavior up to scale 
£* are very disordered so that a weak, elastic-like distortion of such modes is diffi- 
cult to detect. Nevertheless, it has been discovered that the cross-over length is easily 
discernable by eye in the response to a point force or to inflation of a local particle 
dEllenbroek et al., 2006] [EUenbroek et al., 2009a|, as Fig. |03] illustrates. A detailed 



analysis of these data has indeed shown that the response is governed by a crossover 
scale that grows as l/AZ, in accord with (|0.6|) . 



0.2.4 Microscopic criterion for stability under compression 

So far our analysis of the density of states has neglected the effect of pressure on the 
vibrational spectrum. Mundane observations such as the buckling of a straw pushed 
at its tips show that compression affects vibrational modes, and can even make them 
unstable. In general, this effect is important for thin objects, but is irrelevant at small 
strain for bulk solids where plane waves dominate the vibrational spectrum. However 
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f=-k5 



Fig. 0.6 Consider for concreteness a harmonic contact between two particles of length s, of 
stiffness k, carrying a force / = —kS where 5 is the contact elongation. If particles moves 
transversally to the contact by an amplitude x as shown with vertical arrows, the contact 
length increases by an amount proportional to /s following Pythagoras theorem. The work 
produced by the contact force is then proportional to x'^kS/s, corresponding to a stiffness 
fci ~ kS/s = —ke, where e = —5/s is the contact strain, chosen to be positive for compressed 
contacts. 

compression plays an important role if floppy modes are present in the solid, as noticed 
by Alexander in the context of gels ( [Alexander, 1998 ). As we now show, this is also 



true for the anomalous modes introduced above ( Wyart et ai, 2005 



The analysis is based on the well-known observation, see Fig. (|0.6p . that if k is the 
stiffness associated with the longitudinal relative displacements of two particles in con- 
tact, then for transverse displacements the stiffness is ki oc —ke, where e is the contact 
strain. For a pure plane wave propagating in an amorphous solid, the transverse and 
longitudinal components of the relative displacements in the contacts are of the same 
order. The correction induced on the mode energy by the presence of a finite contact 
force or strain is thus, in relative terms, of order of the ratio of the two stiffnesses: 
— e. Therefore the effect of compression on plane waves is negligible at small strain. 
For anomalous modes the situation is different, because they are built by deforming 
floppy modes that have no longitudinal relative displacements, but only transverse 
ones. Once floppy modes are deformed to generate anomalous modes, they gain a lon- 
gitudinal component of order AZ as follows from the variational argument, whereas 
the transverse component remains of order one ( Wyart et ai, 2005 &[ | Wyart, 2005p . 



Relative corrections in the anomalous mode energies is thus of order —e/AZ'^. When 
these relative corrections reach a constant of order one, anomalous modes become un- 
stable, and the system yields. For particles near jamming, one has e oc At/), leading to 
the condition for stability ( [Wyart et ai, 20056 ): 



AZ > A(/)^/2 for all subsystem of size L > I* (0.8) 



where a numerical pre-factor is omitted. Inequality (|0.8p extends Maxwell criterion to 
the case of finite compression and is also not a local criterion: it must be satisfied 
on all subsystems of size larger than /*. On smaller scales, fluctuations of coordina- 
tion violating Eq. (|0.8|) are permitted, as stability can be insured by the boundaries. 
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Configurations where the bound ()0.8p is saturated are marginally stable and anoma- 
lous modes exist down to zero frequency. Such packings should exhibit the scaling 
AZ ~ A0^/^. Configurations obtained by decompression in the vicinity of (j)c ap- 



pear to lie very close to saturation ([Wyart et ai, 20056|. The scaling of Ea. (|0.8p was 



proposed independently by Head (Head, 2005) following a mean- field analysis; there, 



however, the threshold coordination found was half its correct value. 

An argument for why one might expect the system to be marginally stable is 
provided in section 10.6.11 where it will be argued that the realization of the bound 
(|0.8p critically affects the dynamics. The criterion (|0.8|) applies to spheres, but not 
to ellipses where rotational degrees of freedom matter. It can be shown that in the 
later case, compression can have a stabilizing effect if the ellipses make in average 
more contacts where their surface curvature is low. In that case, even an infinitesimal 
pressure can stabilize zero-modes, so that stable packings can be generated although 
the Maxwell bound is not fulfilled, as shown in section ID. 3. II The same is true for gels 



of cross-linked polymers (Alexander, 1998). 



0.2.5 Behavior or elastic constants near jamming 

The fact that a packing of frictionless spheres or discs is a marginally-connected solid, 
also has its effect on the clastic constants. The behavior of the bulk compression 
modulus B and shear modulus G always depends on the specific force law, and hence 
is not universal: for distributions of bonds which are peaked around some average non- 
zero value, both of them are proportional to the average bond strength fc, hence in 
discussing the effects due to jamming it is useful to divide out this common effect. The 
individual bond strength fc^ are essentially the second derivatives of the interaction 
potential evaluated at each contact — for the one-sided harmonic springs these are 
all the same, but for the Hertzian ~ 6^-'^ force law one finds fc.y Sjj^ ^ flj^^ 
that the average bond strength scales with the pressure p as p^^^. 

An example of data for the elastic moduli is shown in Fig. 10.71 which summarizes 
data for B and G for packings of two-dimensional discs interacting with Hertzian forces 
( jEllenbroek et ai, 2006| Ellenbroek et ai, 2009 a| ) . The bulk modulus B is clearly seen 



to scale as p^^^ (O'Hern et ai, 20031, which as we argued above is the scaling of 



the average bond strength for this force law. However, G decreases much faster, 
as p^/^ ( jO'Hern et ai, 2003| ): packings close to jamming are much easier to distort 
by shear than by compression. Results in both two and three dimensions for har- 
monic QDurian, 19"95{|0'Hern et al, 2003| and Hertzian potentials ( jO'Hern et ai, 2003| ) 
have been determined. From these results, it is generally seen that G/B and AZ scale 
in the same way with compression. The data in Fig. 10.71 are consistent with this, as 
for a Hertzian interaction, G/B ~ p^/^ and AZ ^ |A(/)|°'^ p^^^. The result that 
the shear modulus generally scales as G kAZ, can be understood ( j Wyart, 2005D : 
independent of the force law, the ratio G/B should scale as G/B AZ. Empirical 
support for this behavior in grains is reviewed in ( jAgnolin and Roux, 2007| . 

It is interesting to note that at first sight the fact that G/B goes to zero at the jam- 
ming point as AZ, may seem to be the anomalous behavior; however, from the point of 
view of rigidity percolation ( jJacobs and Thorpe, 1995] ) , B/k behaves anomalously for 
packings, as in rigidity percolation both G/k and B/k vanish as AZ, upon approaching 
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Fig. 0.7 (a) As the jamming point is approaclied, in tliis case by lowering tlie pressure, tlie 
siiear modulus G becomes much smaller than the bulk compression modulus B ( [Durian, 1995[ 
|0'Hern et al, 2003| . These results are for packings with Hertzian forces, for which the bond 
strength of individual forces scales as p^^^. Due to this weakening of the bonds, the bulk 
modulus B has an overall p^^^ scaling, while the shear modulus G has an overall p^^^ scal- 
ing, (b) The ratio G/B is linear in Az ( [Durian, 1995| |0'Hern et al., 2003[ l; this scaling is 
independent of the force law. Squares refer to data obtained from response to a local point 
force, diamonds to data obtained from a global deformation. After ( |Ellenbroek et al, 2006| 
[EUenbroek et al, 2009"a| ). 

the isostatic point so that G/B remains constant at the percolation threshold. The dif- 
ference between the jamming transition and generic rigidity percolation is that in the 
former case configurations must be such that all the contact forces must be repulsive. 
This imposes a geometric constraint on the network that is not enforced in standard 
rigidity percolation models ( [Wyart, 2005^ . We refer to ( [EUenbroek et al., 2009 &[ ) for 
further details. 

To compute static properties such as G/B theoretically, one may use effective 
medium theory. This approach has been used to study a square lattice with ran- 



domly placed next-nearest-neiglibor springs (Mao et al., 2010) and random isotropic 
off-lattice spring networks ( [Garboczi and Thorpe, 1985 Wyart, 2010[ ). In the isotropic 
case (rigidity percolation), G/B is independent of AZ at the transition, as found in 
Ref. ( [EUenbroek et al., 2009 6[ ), while in the randomly decorated square lattice G/B ^ 
AZ^. Both results differ from the scaling observed for jammed packings, G/B ^ AZ. 

An important topic in granular research has been the issue of the coarse graining 
scale needed for elastic behavior to emerge in granular packings ( [Goldhirsch and Goldenberg, 2002{ 
[Goldenberg and Goldhirsch, 2008] ). In line with the arguments given above, it has been 
found that near the jamming point, one needs to coarse grain over the scale i* to see 
elastic response emerge (EUenbroek et at, 2006 EUenbroek et al., 2009a|: as one ap- 



proaches the jamming point, one has to coarse grain over larger and larger lengths, 
diverging as £* ~ 1/AZ. That elasticity emerges in granular media on long enough 
scales is also implicitly confirmed by the data of Fig. 10.71 where squares give the 
values of the elastic constants obtained from local point response calculations, and 
diamonds those obtained from global deformations. The two sets of data agree very 



weU (EUenbroek et al., 2006 EUenbroek et al., 2009a). 
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0.2.6 DifFusivity, quasilocalized modes and anharmonicity 



DifFusivity. The excess low-frequency vibrational modes associated with the isostatic 
jamming transition have important consequences for the nature of energy transport 
in the system. Their contribution to energy transport can be quantified by the energy 
diffusivity, defined as follows ( [Allen and Feldman, 1989a| [Allen and Feldman, 1989&I . 
Consider a wavepacket narrowly peaked at frequency uj and localized at r. This 
wavepacket spreads out in time, such that the square of the width of the wavepacket 
at time t, divided by t, is independent of t at large t; this constant defines the dif- 
fusivity, d{uj). The diffusivity was calculated using the Kubo approach within the 
harmonic approximation for repulsive sphere packings as a function of compression in 
( [Xu et ai, 2009[ [Vitelli et ai, 2010| ). At sufficiently low frequencies, the modes should 
be plane-waves and the diffusivity should follow Rayleigh scattering behavior so that 
d{tj) ~ a;-". Refs. ( [Xu et al, 2009] [Vitelh et ai, 2010[ ) show that the diffusivity is 
nearly flat for frequencies from w* up to the Debye frequency, where it drops to zero. 
Thus, there is a crossover from weakly-scattered plane waves to strong-scattering, dif- 
fusive behavior at w*; this can be understood from the heterogeneous nature of the 
anomalous modes ( [Wyart, 2005[ [Vitelli et al., 2010| , or by assuming the continuity 
of the length scale characterizing plane waves and anomalous modes at the cross-over 
frequency, in conjunction with the scaling laws for the shear and bulk moduli discussed 
in section [OTSi ( [Vitelh et ai, 2010[ [Xu et ai, 2009 1 . 

One consequence of the link between the transport crossover and the excess modes ( [Xu et ai, 2009[ 
[Vitelli et ai, 20l"0| ) is that the diffusivity plateau extends all the way down to the low- 
est frequencies studied as the jamming transition is approached and uj* — > 0, as shown 
in Fig. lO.Sf a) . This suggests that the origin of the diffusivity plateau can be traced to 
properties of the jamming transition. 

Dynamical vibrational properties such as the energy diffusivity can be computed 
theoretically using dynamical effective medium theory, a re-summation scheme of a 
perturbation expansion in the disorder also known as the Coherent Potential Approx- 
imation (CPA) ( [Soven, 1969[ ). This approach has been used on the square lattice 
with randomly placed next- nearest- neighbor springs ( [Mao et ai, 2010 1 and isotropic 
disordered spring networks near the isostatic limit ( [Wyart, 2010 1. Both calculations 
yield a crossover to a fiat density of states above uj* ^ AZ in good agreement with 
results for jammed packings. In the isotropic case, Rayleigh scattering is found for 
u) < ui* with a scattering length £s ~ AZ^ /uj'^. This calculation predicts a violation 
of the loffe-Regel criterion (according to which the cross-over to strong scattering oc- 
curs when the scattering length becomes of order of the wavelength) at the crossover 
frequency uj* since the scattering length ig ~ 1/AZ is much larger than the wave- 
length A ^ 1/^/AZ there. Above uj*, the anomalous modes are characterized by a 
displacement-displacement correlation length that scales as l/^/uj, which scales as 
I/VAz at UJ*, and a frequency-dependent speed of sound, c{uj) ^ ^fuj. This leads to a 
plateau in the diffusivity consistent with the results of Fig. (j0.8p . This analysis predicts 
a drop in the speed of sound near uj* , a subtle effect recently observed numerically 
at the Boson peak frequency in model glasses ( Monaco and Mossa, 2009] ) . In the case 
where the isostatic state is a square lattice ( [Mao et al., 2010 ) , the loffe-Regel criterion 
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Fig. 0.8 (a) Difltusivity, d{ijj), for a system of spheres interacting via liarmonic repulsions 
very close to the jamming transition. d(tj) is essentially flat down to zero frequency. From 
Xu et al. ( |Xu et al., 2009p . (b) Energy barrier Vma^ioj) along each mode, before falling into 
another basin of attraction for a system at = 0.1. The energy barriers along each mode 
direction provide an upper bound on the energy barriers of the system in a given frequency 
range, since it is possible that lower energy barriers could be encountered in directions in 
phase space corresponding to linear combinations of the modes. From ( |Xu et al, 2010| . 

is satisfied at the crossover frequency ur* and is is of order the lattice spacing. Numer- 
ical results for jammed packings are consistent virith is ^ l/V Az at the crossover. 

Quasilocalizcd modes. The nature of the excess low frequency modes shows a sys- 
tematic variation with frequency. As the frequency is lowered towards co*, numeri- 
cal studies of three-dimensional sphere packings show that the modes become pro- 
gressively more heterogeneous with a lower-than-average mode coordination number, 



and develop high displacement amplitudes in small regions of space (Xu et al., 2009). 
These quasi-localized modes shift downwards with decreasing compression. These low- 
frequency, quasi-localized modes have also been observed in two-dimensional packings 
of soft disks ( [Zeravcic et ai, 2009[ ), but here the quasi- localized modes are mixed in 
with plane- wave-like modes over the same low-frequency range. 

In three dimensions, the frequency of the quasi-localized modes decreases with com- 
pression, but the degree of quasi-localization does not appear to depend on compres- 
sion. This suggests that as the system is decompressed towards the unjamming transi- 
tion, quasi-localized modes always exist but shift to lower and lower frequency. Quasi- 



localized modes have been observed in other disordered models (Biswas et ai, 1988 



Fabian and Allen, 1997| [Taraskin and Elliott, 1999| [Schober and Ruocco, 2004[ ) and 

in experiments on glassy polymers ( Buchenau et al, 1992[ Buchenau et al, 1996| [Vainer et al, 2006D . 

In the soft-potential model, for example, such modes have been shown to lead to strong 

scattering of vibrations at higher frequencies, so that their frequency corresponds to the 

loffe-Rcgel crossover frequency in the transport properties ( [Schober and Oligschleger, 1996[ ), 

while in model silica, these modes have been identified near the boson peak fre- 



quency (Taraskin and Elliott, 19991 



Numerical studies on hard spheres ( [Brito and Wyart, 2007[[Brito and Wyart, 2009[ ) 
have shown that irreversible particle rearrangements at nonzero temperature occur 
along the low- frequency modes, observed to be extended in the systems of limited size 
studied, implying that the lowest activation barriers are in those directions of phase 
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space (Sec. l0.6.l|) . This view has been explored in model glasses dWidmer-Cooper et ai, 2009 aD , 
which have also shown that the structural rearrangements occurring at finite tempera- 
tures were along low-frequency normal modes, which in that case were quasi-localized. 

Fig. IO.8f b) shows that for sphere packings, the low- frequency quasi-localized modes 
have the lowest energy barriers and that the energy barriers decreases systematically 
with frequency and participation ratio. This result, together with previous observa- 
tions dBrito and Wyart, 2007] |Brito and Wyart, 2009] [Widmer-Cooper et ai, 2009aD , 
should have important consequences for the dynamics and instabilities that occur as 
the temperature is raised from zero, the system is compressed, or a shear stress is 
applied. One would expect these modes to shift downwards to zero frequency under 



the last two circumstances ( Maloney and Lemaitre, 2006 ) and to give rise to highly 



heterogeneous deformations in the sample. These could be the underlying mecha- 
nism for the localized deformations postulated in the picture of shear transformation 
zones dFalk and Langer, 1998| [Langer, 2008] ). 

As the system is decompressed towards the unjamming transition, the barriers 
between nearby configurations shrink to zero as the amount of inter-particle overlap 
decreases to zero. Thus, in this limit the modes become progressively more anharmonic. 
This implies that smaller excitation energies or temperatures should suffice to force 
the system into new energy minima. Thus, the range of temperatures and stresses 
over which harmonic theory applies should shrink to zero and anharmonic effects 
should dominate as one approaches the transition. This effect has been computed in 
networks of springs below or above the isostatic threshold ( |Wyart et ai, 2008] ) , where 
non-linearities dominate the response to shear when the strain 7 is larger than some 
7*, with 7* - |Az|. 

0.2.7 Structural features of the jamming transition 



At point J, the pair distribution function has many singular properties ( O'Hern et ai, 2003 
[Silbert et ai, 2006| . In particular, it has been shown that the first peak in g{r) is a 
(5-function at the nearest neighbor distance, a. The area under this (5-function is just 
the isostatic coordination number: .Z = 6 in three dimensions. On the high side of 
this i5-function, g{r) has a power-law decay: g{r) cx (r — a)~^-^ . It has been proposed 
that this power-law is the vestige of the marginal stability of the configurations visited 
before reaching (j)c d Wyart, 2005D . Also, at the jamming transition, the second peak of 
g{r) is split into what appear to be two singular sub-peaks one at r = -s/Sc and the 
other at r = 2a. There is a divergent slope just below each subpeak and a step-function 



cutoff on its high side( Silbert et ai, 2006). 



The pair correlation function near the zero-temperature jamming transition has 
been studied in two experiments. In experiments by Abate and Durian dAbate and Durian, 2006} 
[Keys et ai, 2007D , a layer of ball bearings were placed on a wire mesh and excited by 
the upflow of gas through the mesh. The turbulent flow of the gas leads to interactions 
between the balls and to random motion of the balls. Thus, the effective interaction 
potential between the balls can be characterized by a hard core determined by the ball 
diameter and a longer-ranged repulsion due to the gas. In such systems, the kinetic 
energy of the balls decreases monotonically with increasing density of the balls, such 
that it reaches zero when the balls are packed to the density of the jamming transition. 
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Fig. 0.9 Pair-correlation function g{r) for tlie large particles at different packing fractions (p 
for a bidisperse system of colloids in two dimensions. The inset shows gi, the height of the first 
peak of g{r), versus packing fraction. The experimental data shows that gi is non-monotonic 
with a maximum at a packing fraction slightly above (f>c, the value at the T = jamming 
transition. The peak is a vestige of that T = transition where the first peak of g{r) is a 
5-function. Figure reproduced from ( |Zhang et al., 2009] ). 

0c — 0.84, measured in terms of the hard-eore diameter. The structure of the system, 
as characterized by the pair correlation function, g{r), was studied along this trajec- 
tory to the zero-temperature jamming transition of the hard cores (See the chapter 
in this book by Dauchot, Durian and van Hecke for a description of the dynamical 
properties of this system). On approaching (f>c, the structure shows an increase of the 
height of the first peak, gi. This is consistent with numerical results on soft spheres, 
which find a divergence in gi at the zero-temperature jamming transition. There is 
also a local maximum in gi at (/> sa 0.74, above which no changes in the coordina- 
tion and geometrical features of Voronoi cells are observed. This secondary maximum 
at nonzero kinetic energy due to airfiow- induced interactions, may correspond to the 
finite-temperature vestige ( [Zhang et ai, 2009| of the divergence in the height of the 
first peak of g(r) at the zero-temperature transition discussed below. 

Experiments on a two-dimensional bidisperse system of NIPA particles by Zhang, 
et al. ( [Zhang et ai, 2009| , measured the pair correlation function as a function of den- 
sity at essentially fixed, nonzero temperature. The system exhibits a maximum in the 
height of the first peak of the pair correlation, gi — see Fig. 10.91 This maximum is 
interpreted as a thermal vestige of the zero-temperature jamming transition. Corrobo- 
rating simulations show (Zhang et al., 20091 that the divergence in gi ~ 1/(0 — (/)c) for 
soft repulsive spheres at the zero-temperature jamming transition is softened to a max- 
imum at nonzero temperature. The maximum in gi decreases in height and shifts to 
higher density with increasing temperature. The simulations also show that the maxi- 
mum in gi is a structural feature that is apparent as a function of increasing density at 
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fixed temperature, but not as a function of temperature at fixed density or pressure. 
Recent simulations on tlie same model by Jacquin, et al. ( [Jacquin and Berthier, 2010[ ) 
show that the maximum in 171 appears even in systems at higher temperature or lower 
pressure that are in thermal equilibrium and can be understood from liquid state 
theory. 

Note that the divergence in gi corresponds to a vanishing of the overlap distance 
between neighboring particles. Thus, tracking gi is equivalent to tracking the overlap 
distance. The evolution of the overlap distance with temperature is not correlated with 
the existence of a dynamical glass transition. 

0.2.8 Singular length scales at Point J 

Several different singular length scales have been identified near the jamming transi- 
tion. The cutting length, £* , discussed in Sec. 1.2.3, has been observed by studying 
the fluctuations in the change of the force on a contact as a function of the distance 
of the contact from a small local perturbation (see Fig. 1.5). This length scale di- 
verges as (0 — (/)c)~^^^, independent of the potential as long as it is repulsive and 
finite in range. In addition, this length scale is expected to correspond to the wave- 
length of longitudinal, weakly-scattered plane waves just below the crossover frequency 
ui* ( Silbert et ai, 2005] ), although it has not been seen numerically, due to the small 



system sizes that can be studied. 

A second diverging length scale has been identified from the wavelength of trans- 
verse, weakly-scattered plane- waves just below the crossover frequency ui* ( [Silbert et ai, 2005{ 
|Vitelli et al., 2010| ). This length scale diverges as (0— 0c)~^^^- Finally, this length scale 
emerges within the effective medium approximation as the decay length for spatial cor- 
relations in modes at the frequency ui uj*+ ( [Wyart, 2010[ ). 

As noted in Sec. 1.2.5, there is a vanishing length scale corresponding to the 
overlap distance S between neighboring particles. This distance, which measures the 
Icfthand width of the first peak of the pair-correlation function, vanishes as {(f> — 
dO'Hern et al, 2003[ [SUbert et al, 200E^ . 



Finally, there are other diverging length scales whose origins are not understood. 
The shift of the position of the jamming transition, 0c, with the linear size of the 
system, L, yields a length scale that apparently diverges as \4> — 0c l""'"^ in both 
2 and 3 dimensions ( [O'Hern et al., 2003[ ). This scaling shows up in simulations in 
which a hard disk is pushed through a packing below (jjc', the transverse distance 
over which the packing adjusts as the disk passes by diverges with a power-law of 
0.7 ( [Drocco et al., 2005[ ). Finally, this power- law has been observed for correlations of 
the transverse velocity on athermal, slowly-sheared sphere packings near the jamming 
transition ( [Olsson and Teitel, 2007 ). However, this length scale is only observed for 



certain models of the dynamics ( Tighe et al., 20101 ) 



0.3 Extensions of the results for frictionless spheres 

0.3.1 Anisotropic particles: ellipsoids 

The counting of degrees of freedom that underlies the above analysis of the soft modes 
near the jamming point, relies on having perfect spherical and frictionless particles. 
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We briefly review here what happens to the jamming scenario if the particles are 
frictionless but non-spherical. 

The model for which this question has been studied in detail is that of two- 
dimensional frictionless ellipses in two dimensions ([Mailman et ai, 2009 ) or ellipsoids 
with one axis of symmetry in three dimensions ( [Zeravcic et ai, 2009 1. In the latter 
case, each ellipsoid is essentially characterized by 5 nontrivial degrees of freedom, so 
a naive counting of the degrees of freedom shows that the isotropic values for such 
ellipsoids is ZfJ^'P-'d ^ 2 • 5 = 10. 

Clearly, while the isotropic values jumps depending on which degrees of freedom 
are, and which ones are not, included in the counting, one would expect the physical 
behavior of such packings to evolve continuously, when the a-sphericity of the particles 
is turned on — what is happening? 

Let us take the c axis of our ellipsoids along their symmetry axis; the other axes 
are then a and b ~ a. We define the ellipticity as the aspect ratio s = c/a, so 
prolate ellipsoids (like cigars) correspond to £ > 1 while oblate ones (like M&M's 
(Donev et al., 2004)) to e < 1. Fig. lO.lOl b) recovers the finding of ponev et ai, 2004[ 



[Donev et ai, 2007| Sacanna et ai, 2007| [Wouterse et ai, 2007[ ) that as the ellipticity 
is tuned away from the spherical case £ = 1, the average contact number indeed 
increases continuously from the spherical value Ziso = 6, as expected. 

How to resolve this apparent paradox between the jump in the counting and the 
continuity of the physical system? How come one can reach values of the average 
contact number below the isostatic value ^°'^'p*'°"^ = iq? Xo answer these questions, 
it is good to realize that when we reached the conclusion that the isostatic value 
^sphere ^ 2d = 6, we left out the rotational degrees of freedom of each sphere, as 
these arc trivial zero-frequency modes (three per particle) for a system of frictionless 
spheres. In order to understand the continuity, it is better to include these rotational 
degrees of freedom from the start, even for the spherical case, and see how they evolve 
when the spheres are deformed into ellipsoids. Since in the work ( [Zeravcic et ai, 2009[ ) 
that we will summarize below, the ellipsoids still have one symmetry axis, there are 
only two nontrivial angular degrees of freedom associated with the anisotropy of these 
particles. The relevant isostatic value for these ellipsoids of revolutation is therefore 
■^is'o^^°"^ = 2 • 5 = 10, as already asserted above. 

When analyzed more carefully, the counting arguments summarized in section lO.2.21 
imply that whenever Z < Ziso, there are Nc{Z{so — Z)I2 directions in the phase space 
of the relevant coordinates of the particles along which the packing coordinates can 
be changed without changing the interaction energy of the particles. This means that 
each stable configuration must have precisely this number of zero-frequency modes. 
That this is precisely what happens in the ellipsoid packings, is illustrated by the grey 
crosses in Fig. lO.lOT b): numerically it is found that in the calculation of the density 
of states of ellipsoid packings one finds precisely the right number of zero modes as 
expected according to this argument. 

Fig. lO-lOf a) illustrates how, as Z increases with Se ~ e — 1 oi the ellipsoids, the 
number of zero modes per particle with contacts goes down. We thus have more 
and more nonzero modes as Z increases. What is their nature? Fig. lO.lOf a) shows 
that for \5e\ < 0.17, the density of states exhibits two clearly separated bands 
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Fig. 0.10 (a) Illustration of ho-w the number of different modes per ellipsoid (excluding 
rattlers) in the packings at jamming varies as a function of the average contact number Z. 
As the ellipticity is increased, the average contact number increases. For Z <« 9, there 
are two well-defined bands: a rotational band with {Z — 6)/2 modes per particle and a 
translational band with three modes per interacting particle as is the case for spheres (see 
Fig. lO.lil for details). Upon increasing Z, the number of zero modes decreases as zero-modes 
are converted into finite-frequency rotational modes. Above Z ~ 9, there is only one band, 
(b) Z as a function of the ellipticity Se = e—1 and distance from jamming 5(f) for 216-particle 
packings. The sharp decrease around the spherical case 5£ = is consistent with earlier 
results ( |Donev et al., 2004[ |Sacanna et al., 2007[ |Wouterse et al., 2007[ l for hard ellipsoids 
and spherically capped rods. The log-log plot of SZ versus \5e\ in the inset shows that the 
rise of Z at jamming is consistent with a SZ ~ -v/l&T scaling. The crosses in the main plot 
for small values of Se show that twice the measured number of zero-frequency eigenmodes 
per particle plus Z — 6 add up precisely to 4, in accord with panel (a). From Zeravcic et al. 
( |Zeravcic et al, 2009| . 

( [Zeravcic et al, 2009p . The modes in the upper band of the density of states turn 
out to be predominantly translational — they are very similar to those found in the 
sperical packings, and turn out to be characterized by an onset frequency w* that 
scales as AZ = Z — 6 ^ \Se\^^^: even though contacts at jamming are created by 
deforming the particles and not by compression, the mechanism that determines the 
onset of anomalous translational modes at a frequency oj* is the same! Furthermore, 
the lo"wer frequency band that exisits for small ellipticity consists of predominantly 
rotational modes; it is characterized by an upper frequency that scales linearly in |fc| 
( [Zeravcic et al., 2009D . For excentricities larger than about 0.17, the t-wo bands merge 
into one mixed band, as Fig. lO.lOf b) illustrates. 

The scenario that emerges for the number of modes in the various bands is sum- 
marized in Fig. 10. Uf a) — the upshot is that the jamming scenario as formulated for 
frictionless spheres, as exemplified by the sharp onset of anomalous modes, "withstands 
the change in shape of the particles in the system. 

0.3.2 Packings of frictional spheres 

Another relevant extension is to include static friction in the force la"w. The usual 
friction la"w is to take the Mohr-Coulomb la"w pohnson, 1985[ ), in "which the particles 
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Fig. 0.11 (a) The density of states for sliglitly oblate ellipsoids, Se = —0.04, for pack- 
ings close to jamming and at two compressions. The existence of two bands separated by a 
gap as well as a gap at zero frequency is clearly visible, (b) For larger ellipticities the two 
bands merge, as is illustrated here for Se — —0.33. The inset shows the data for Se = 1 
where Z « Zc = 10 at jamming (see panel a). In accord with this, the gap near zero fre- 
quency increases with increasing compression (and therefore increasing Z). From Zeravcic et 
al. ( | Zeravcic et al, 200"9| . 

can experience a tangential force ft which obeys ft < ^J■fn, where /„ is the normal force 
at the contact considered, and ^ the Coulomb constant. In the Hertz force-law model, 
the tangential elastic force builds up continuously, in a well-described way, when a 
contact gets loaded with a tangential force ( [Johnson, 1985 ). This buildup continues 
up to fxfn, beyond that ft remains constant and fixed at the value nfn- In the presence 
of static friction, therefore, the properties of a packing are very history dependent, as 
the tangential forces are not only dependent on the configuration, but also on how 
the configuration was made — see e.g. ([Silbert et al., 2002{ [Zhang and Makse, 2005 



Somfai et al., 2005 Somfai et al, 2007). 



If we return to spherical particles, but include the Mohr- Coulomb friction law, 
again the counting of the number of degrees of freedom is different: the coordinates 
of the centers of the particles are still the degrees of freedom with which one can 
satisfy the "just touching conditions" (j0.2p upon approaching the jamming point. 
On the other hand, the counting of the degrees of freedom in the force and torque 
balance conditions is different — e.g., in two dimensions the frictional forces add one 
additional tangential frictional force component for each contact, but torque balance 
only adds one extra constraint per particle. A straightforward generalization (see e.g. 
( [Shundyak et al., 2007P ) shows that in general with frictional forces stable packings 
can be obtained provided that d + 1 < Z < 2d. Unlike for the frictionless case, 
where this argument dictates Z completely in the limit that the jamming point is 
approached, the inclusion of a static frictional force allows a range of possible values 
of Z. In simulations of frictional systems, it is indeed indeed found that packings 
at the jamming threshold with different values of Z and density (j) can be made, by 
varying in particular the quench rate; moreover, there is a clear trend for the average 
contact number Z to approach the lower lower limit d + 1 — the isostatic value in the 
presence of friction — in the limit ^ — >• oo ( [Shundyak et al., 2007[ [Song et al., 2008P . 
In line with this, one finds that if one studies the linear vibration-rotation modes of 
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Fig. 0.12 Data from Menkes et al. ( |Henkes et al, 20106P on frictional packings, (a) In gently 
prepared frictional packings, there a finite fraction of the contacts is found to lie at the 
mobilization threshold (are "fully mobilized"), i.e., has ft — iifn- In this diagram a large 
number of packing configurations are plotted in a diagram of Z and rim, the average number of 
fully mobilized contacts per particle. For lower pressure, the dots in this diagram representing 
the packings approach the generalized isostaticity line, defined by (|0.9p . where packings are 
marginal, (b) Density of states for the lowest pressure data for a large range of values of 
the friction ^. In all cases the density of states is strongly enhanced at low frequencies. The 
frequency plotted here is actually the scaled frequency uj mentioned in section [0.2.31 (c) The 
cross-over frequency uj* scales as /S.Z = Z — Zl^o- 

these packings, the density of states only develops a plateau in this large friction limit 
dSomfai et al, 2007D . 

However, it was later discovered that the situation is more subtle. If packings are 
prepared gently, then one finds that a finite fraction of the contacts in the simulations 
are at the mobility or Coulomb threshold, i.e., have ft — fifn ( |Shundyak et al., 2007 



Henkes et al., 20106[ ) . If such a contact is pushed in the direction of increasing friction 



force, they slip! In other words, a finite fraction of the tangential forces is not arbitrary, 
but constrained to be at the Coulomb threshold. This of course changes the counting 
of the number of degrees of freedom. Indeed, if one introduces rim, the number of 
contacts per particle which is found at the mobility threshold, then one may introduce 
a generalized criterion ( [Shundyak et al., 2007 Henkes et al., 20106 1 



9n 

Z>(d+1) + -^^Z™. (0.9) 

for the stability of packings. In this scenario, packings which are close to the line 
Z = Zj™ = d + 1 + 2nm/d are again marginal, i.e., develop many low-lying excita- 
tions because the number of contacts that slip under tangential motion is just large 
enough to make almost zero-energy deformations possible. The idea that gently pre- 
pared packings might end up being marginal packings in this way, traces back to an 
earlier suggestion by Bouchaud ( [Bouchaud, 2004[ ) . 

With this generalization, packings with static friction that are gently prepared, do 
again follow in essence the scenario described above. As illustrated in Fig. l0.12T a). as 
the pressure on packings is reduced, the points representing these packings in a Z — 
diagram approach the generalized isostaticity line ( Henkes et al., 20106[ ). Moreover, if 



the contacts at the mobility threshold are allowed to slip in analyzing the rotation- 
vibration dynamics, the density of states of packings close to the jamming threshold 
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is strongly enhanced at low frequencies for all values of the friction coefficient /x — 
see Fig. I0.12f b'). Moreover, as the pressure is increased to tune the packings away 
from the generalized isostaticity line, the crossover frequency w* in the density of 
states scales linearly with AZ = Z — Z;^, i.e., linearly with the distance from the 
generalized isostaticity line. This is the first confirmation of the robustness of the 
concepts reviewed above to frictional particles. 

There are still several unresolved subtleties associated with the somewhat singular 
nature of the Coulomb threshold condition, but we refer for a discussion of these to 



(Henkes et ai, 201061. Another interesting recent development (Henkes et ai, 2010al 
is that precursors to avalanches in simulations of tilting of piles of frictional spheres 
are found in the correlation properties of balls with low contact number. 

0.3.3 Finite shear forces: nontrivial rheology in the sheared bubble 
model 

So far we have only discussed static packings and their static and dynamic linear re- 
sponse. An important development of the last two years is that simulations of bubble 
models under finite shear stress or at a constant shear rate are probing the jam- 
ming phase diagram along the stress direction. The bubble model ( [Durian, 1995[ 
|Durian, 1997[ ) is essentially a soft-sphere model, enriched with viscous friction, i.e., 
with dynamical terms proportional to the difference between the velocity of a bubble 
and that of each one it makes contact with. A particular attractive feature of such 
bubble model simulations with dynamic friction is the fact that the relevant reference 
frame for the static limit is the frictionlcss sphere or disc model, about which so much 
is known. 

Olsson and Teitel ([Ol sson and Teitel, 2007[ ) were the first to demonstrate that the 
rheological data for shear flows in such systems show very good data collapse close 
to the jamming point of the frictionless case. Fig. 10.131 shows a result from their 
simulations. In this particular plot, the effective viscosity 77 = ct/7 is plotted as a 
function of the shear stress a imposed on the system; here 7 is the shear rate. By 
scaling both quantities properly with the distance At/) from the jamming density, there 
is very good data collapse both above and below the critical density. Moreover, at the 
jamming density, and close to the jamming density for large enough shear rates 7, the 
data is consistent with a scaling <t ^ 7" "'. The excitement in the field is due to the 
fact that this is one of the first examples of how a simple microscopic model leads to 
non-trivial rheological behavior at mesoscopic or macroscopic scales, and the fact that 
this is intimately related to the jamming scenario. 

The more customary way to plot rheological data is to plot the stress cr as a function 
of the imposed shear rate 7; this is how data like that of Fig. lOTBl are usually presented. 
At the time of writing this review, there is a lot of debate about the precise exponents 
(Olsson and TeitcQlWllHatano, 20086||Hatano, 2008a|[H"atano, 2008c||Olsson and Teitel, 20091 



ighe et ai, 20lO[|Haxton and et ai, 2009 1, the existence of a quasistatic limit, the in- 



fluence of the mean- field approximation employed in the simulations of ( Olsson and Teitel, 2007 ) 
the question of whether there is a finite velocity correlation length which diverges at 
the jamming point, or one that is always of order the system size, etc. From our point 
of view, the essential question however is whether, just like in the static case, the expo- 
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Fig. 0.13 Scaling plot of Olsson and Teitel ( |01sson and Teitel, 2007[ ) for simulations of the 
bubble model under shear. Plotted is a the viscosity rj = cr/7, rescaled by the density, as a 
function of the applied shear stress cr, with 7 the shear rate. The scaling collapse illustrates 
that close to the jamming point, nontrivial rheology occurs. Close to the jamming point, we 



ncnts should be rational values which can be understood in terms of simple arguments 
based on our insight in the response of a marginal solid. To date, the answer to this 
question is unknown. 



0.4 Real physical systems 

As stated in the introduction, space limitations do not allow us pay tribute appropri- 
ately to the recent experimental developments that provide new impetus to the field. 
Some of these are reviewed in the chapters by Cipelletti and Weeks on colloids and by 
Dauchot, Durian and van Hecke on grains and foams. We wish to draw attention to a 
few additional points here. 



0.4.1 Granular materials 

Excited granular matter turns out to be an excellent model system to probe dynami- 
cal heterogeneities and the pair-correlation function g{r) — think of the experiments 

Lechenault et ai, 2008 6p or 



with discs shaken on a table (Lechenault et al, 2008a 



the ball bearing experiments ( [Abate and Durian, 2006 



Keys et al., 2007) discussed 



in section [p. 2. 71 Unfortunately, it is extremely difficult to probe the vibrational den- 
sity of states with grains. The scaling of static properties like the contact number is 
also difficult to study experimentally: in experiments one will always have friction, 
whereas most of the clean theoretical scaling predictions are for frictionless spheres 
— as we have seen in section 10.3.21 with friction the behavior depends stongly on 
the way a packing is made, so it is not even clear whether clear scaling behavior 
should be expected in the experiments on photoelastic discs ( [Majmudar et ai, 2007[ ) 
aimed at measuring Z as a function of packing density (However, experiments using 
photoelastic discs have been very powerful in studying force-chain networks and cor- 
relations ( [Majmudar and Behringer, 2005[ ).) A powerful way to probe elasticity is to 
use ultrasound measurements, which, in agreement with what we have argued above, 
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support a small G/B ratio in sand ( [Bonneau et ai, 2008[ [Jacob et ai, 2008"! ) and a 
sharp transition toward strong scattering ( Jia et al, 1999 ) as the frequency increases. 
Another possibility which is being explored to probe the contact statistics and elas- 
tic moduli by using particles which swell when the humidity of the environment is 
changed ( [Cheng, 2010[ ) . Finally, experiments on avalanches are convenient for probing 
the length scales involved in failure and flows ( [Forterre and Pouliquen, 2008[ ). In par- 
ticular the angle at which an avalanche starts depends strongly on the thickness of the 
granular layer. A model based on the stabilization of anomalous modes by the rough 
fixed boundary at the bottom of the granular layer is consistent with the observed 
dependence ( [Wyart, 2009 a[ ). 



0.4.2 Emulsions and foams 

It is important to stress here that some of the observations associated with the jam- 
ming approach — in particular, the scaling AZ ~ \/ Acj) — already appeared in some 
form in the early work on bubbles and the bubble model ( [Durian, 1995[[Durian, 19"97| . 
Moreover, it has become increasingly clear that wet foams are a very good experimen- 
tal model system for jamming studies, since the friction is dynamic, not static. As a 
result, models of frictionless spheres or discs provide a reference frame to study the 



effect of shear. Indeed, recent experiments on sheared foams (Katgert et ai, 2009& 



Katgert et ai, 2009a) hold the promise for testing some of the predictions of the 
sheared bubble model that were discussed in section [CT.3. 31 

Another promising route is to use emulsions or foams. Here too, experiments from 
over a decade ago ( [Mason et ai, 1995\ [Saint- Jalmes and Durian, 1999[ ) already gave 
indications for deviations from the behavior expected for affine deformations. New, 
more-refined experiments, are called for. Emulsions have recently also been used to 



probe the structure of static packings (Clusel et al., 20091 



0.4.3 Colloids 

Soft colloids, such as foams, emulsions and particles made from poly(N-isopropyl acry- 
lamide) microgels (NIPA particles) ( [Saunders and Vincent, 1999} [Pelton, 2000[ ), can 
be used to study densities above the jamming transition. In aqueous solution, NIPA 
colloids swell substantially as the temperature is reduced only slightly. As a result, 
the packing fraction of the sample can be varied over a wide range with only a min- 
imal change of temperature. The systems are thermal in that the particles are small 
enough (micron-sized) to exhibit significant Brownian motion. Thus, such systems are 
well-suited for studying jamming at a fixed non-zero temperature. 



0.5 Connection with Glasses 

In this section, we will describe some of the potential connections between the physics 
that appears at the jamming transition at point J and what happens in real molecular 
glass formers. There are several phenomena that appear ubiquitously in amorphous 
solids and glasses yet are quite different from what is found in ordinary crystals. 
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0.5.1 Properties of glasses to be considered 



Density of low-temperature excitations and the Boson peak: There are many more low- 
frequency excitations in a glass than in a crystal. One of the hallmarks of glasses are the 
so-called anomalous low-temperature excitations. Since 1971 ( Zeller and Pohl, 1971[ ), 
it has been known that the specific heat, cy of glasses varies approximately linearly 
with temperature rather than as as would be predicted by the Debye calculation 
where the low-frequency modes are long- wavelength planewaves( |Ashcroft and Mermin, 1976[ ). 
At a somewhat higher frequency, above a frequency ujbp-, there is a dramatic excess 
of vibrational modes in glasses, known as the Boson peak(Sokolov et ai, 1995). That 
the Boson peak should be a ubiquitous property of amorphous solids is surprising 
since the underlying structure, which would normally determine the phonon spectrum 
in a crystal, is so different in different glasses. The linear specific heat is surpris- 
ing since even disordered systems should have well-defined elastic moduli when aver- 
aged over large enough wavelengths. Therefore, at low enough frequencies, the modes 
should remain planewaves with a Debye spectrum. This linear specific heat has been 
ascribed to the quantum tunneling of groups of atoms between roughly equivalent 
structural configurations ( [Anderson et ai, 19"72{ |Phillips, 1972[ ) although these config- 
urations have not been unambiguously identified. 

Thermal conductivity: In a related vein, the heat transport in glasses is decidedly 
different from what it is in crystals! Anderson, 1981] ). At very low temperatures, the 
thermal conductivity, k ^ rather than k ~ T'^ as in crystals. At higher temperatures 
there is then a plateau region. Above the plateau there is a third regime where k 
increases gently, /t ~ T. This behavior is surprisingly universal between different glasses 
yet there is no generally-accepted explanation of either of the higher-temperature 
regimes. The lowest-temperature region, where k ~ T^, has been modeled as the 
scattering of plane waves off the localized tunneling systems used to explain the low- 
temperature linear term in the specific heat. 

Failure under applied stress: When an amorphous solid is compressed or sheared 
sufficiently, it will start to fail. However, the failure often occurs in the form of rather lo- 



calized rearrangement events rather than via a long-range catastrophic collapse! Falk and Langer, 1998 



[Langer, 2008[ Widmer-Cooper et ai, 2009a|. Since the system goes unstable when the 
frequency of a mode reaches zero, it is possible to identify a rearrangement with the 
lowest- frequency mode just before the rearrangement ( [Maloney and Lemaitre, 2006[ ). 
However, ideally one would like to predict well in advance where failure will occur, 
even in large systems with many plane- wave-like modes at low frequencies. 

Structure: In metallic and colloidal glasses, the structure has often been described 
in terms of sphere packings ( [Finney, 1977| [Cahn, 1980[ [Zallen, 1998[ ). While this ap- 
pears to be quite natural, there are certain features that have been related to the 
glassy structure that are not well understood. For example, the first peak in the pair 
distribution function, g{r) is tall and sharp and there is a split second peak. Moreover, 
as the temperature is lowered in a supercooled liquid, there is no sign of the onset 
of rigidity in the static structure factor, S{q); the first peak in S{q) simply grows 
smoothly in height and decreases in width as the temperature is lowered through the 
glass transition temperature, Tq( [Busse and Nagel, 1981[ ). 
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0.5.2 Relating the Jamming Paradigm to Glasses 

Low-temperature excitations and the Boson peak in covalent glasses 

Various theoretical approaches can reproduce a Boson peak and some of the fea- 
tures of the corresponding normal modes. They are based on disorder and generally 



connect the peak to an elastic instability, as in effective medium models ( Schirmacher et al., 2007 



[Schirmacher, 2006 1, euclidian random matrix theory ( [Grigera et al., 2001 1 or mode 



coupling ( Gotze and Mayr, 2000[ ). However, the structural parameters controlling the 



peak are imposed by hand in these theories (amount of disorder in ( Schirmac her, 2006] ) 



density in ( [Grigera et al., 2001[ ) or structure factor in ( [Gotze and Mayr, 200^ ). What 
aspects of the structure are the most relevant for the peak is unclear. Moreover, dis- 
order cannot be a generic explanation for the Boson peak: silica has one of the largest 
peaks, but so do the corresponding crystals, the cristobalites, as shown in Fig. 10.141 It 
is sometimes argued that silica is a special case. The variational argument of section 
I0.2.3l shows that is not so: generically amorphous and crystalline elastic networks must 
have qualitatively similar low-frequency spectra if their coordination is identical. For 
example, a cubic lattice has a flat spectrum as does a random close packing of elas- 
tic particles. Thus coordination rather than positional disorder matters (but disorder 
in the spring stiffnesses obviously matters and can be incorporated in an improved 
variational argument ( [Wyart, 2005[)). 



The variational argument of (Wyart et al, 2005a) has broad applications, as illus 



trated for silica and other tctrahedral networks following the analysis of ( [Wyart, 2005[ 
[Wyart, 2009 6[ ) . A particularity of silica or germanium oxide is the weak force constant 
associated with the rotation of two adjacent tetrahedra, in comparison with deform- 
ing the tetrahedra themselves. This suggests modeling silica as an assembly of stiff 
tetrahedra connected by flexible joints, the RUM model ( [Tkachenko et al., 2000[ ), as 
sketched in Fig. 10.141 When this approximation is used on configurations obtained at 
various pressures to compute the vibrational spectrum, one finds a remarkable simi- 
larity between silica and particles near jamming as illustrated in Fig. 10.141 The cause 
is identical: at reasonable pressures, silica is made of tetrahedra, and 5-fold defects 
are extremely rare. If joints are flexible, a tctrahedral network is isostatic (there are 
6 degrees of freedom per tetrahedron, and 3 constraints per joint, and twice as many 
joints than tetrahedra) and must have a flat density of states, as observed. For real 
silica, the rotation of the joints has a finite stiffness, which must shift the anomalous 
modes forming the plateau by some frequency scale A. Using ab initio computation 
of the joint-bending force constant, and the molecular weight of silica one estimates 
A = l.AThz. Empirically the density of states of silica indeed has a plateau above 
approximatively 1 Thz, as expected. This is shown in Fig. 10.141 for numerical simu- 
lations. This argument also applies to germanium oxide. For amorphous silicon, the 
bending force constant of the joints is comparable to the other interactions, since they 
correspond to the same covalent bond, and one therefore does not expect a Boson peak 
at low frequency; indeed in this material the Boson peak is nearly non-existent. Ref. 



(Wyart et al., 2008) extends these arguments to other covalent networks. 
Low-temperature excitations and the Boson peak in molecular glasses 
We can make progress in understanding the excess number of low-frequency exci- 
tations in molecular glasses by recalling some of the properties of solids near the jam- 
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Fig. 0.14 Left: Configuration of amorphous silica as modeled by RUM: Si04 tetrahe- 
dra are connected by springs of finite stiffness at tfie joints, whicli can freely rotate 
( |Tkachenko et al, 2000[ ). Center: Density of states D{uj) for amorphous silica configurations 
obtained at various pressure within the RUM approximation. At low-pressure, there are only 
tetrahedra, and D{ijj) is essentially fiat. As the pressure increases 5-fold coordinated Silicon 
atoms become more and more numerous, and D{ijj) erodes at low-frequency due to increased 
coordination, very much like in sphere packings ( [Tkachenko et al, 2004[ |. Right: Computa- 
tions of D{uj) for atomic models of silica. Silica glass, /3-cristobalite and a-cristobalite all 
present a plateau above 1 Thz, the Boson peak frequency ( [Tkachenko et al, 2000[ ). 

ming transition. At point J, the jamming transition for frictionless repulsive spheres, 
the density of states has a plateau all the way down to zero frequency. Upon compres- 
sion, the plateau persists only down to a lower cutoff frequency, a;*( Silbert et ai, 2005[ 



lO'Hern et al, 2003P [Figs. and ES^b)] as explained in ( |Wyart et ai, 2005 aD . Be- 
low that frequency, the normal modes are weakly scattered plane waves while above 
it they are strongly scattered plbert et al., 2009| |Xu et al, 20091 IVitelh et ai, 2010 



[Wyart, 2010[ ) . This jump in D{uj) appears to be the counterpart to the dramatic in- 
crease in the number of vibrational modes that appear at the Boson peak in glasses( [Silbert et ai, 2005[ 
|Xu et al., 2007||Wyart, 2005] [Wyart et al, 2005a[ ) . 

In order to assert this correspondence more generally, we have ascertained that the 
physics dominating point ,J for repulsive spheres is also operative for other systems 
such as those with particles of non-spherical shape or with particles mediated by long- 
range and/or attractive interactions. To include such effects, we have studied several 
different models of glasses that can be analyzed in terms of their proximity to the 
jamming transition. We refer to section lD.3.1l on ellipsoids for a brief discussion of what 
happens when the particles are non-sphcrical — in essence, the results are consistent 
with the jamming scenario based on the analysis of frictionless spheres, though in a 
surprising and nontrivial way. 

We have also studied models that were based on variants of the Lennard- Jones 
interactions between spheres. Such potentials allow attractive as well as long-range in- 
teractions so that the average coordination number can be arbitrarily large. Lennard- 
Joncs potentials are more realistic models of molecular systems than the simple har- 
monic force law described above because the potential decays rapidly with inter- 
particle separation. Particles separated by even a slightly larger distance, interact 
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much more weakly than those closer together. 

We have shown ( jWyart, 2005[ |Xu et al., 2007[ ) that the onset of excess modes in 
these systems derives from the same variational-argument considerations that arise 
at point J for systems with finite-ranged repulsions. It can again be estimated by 
analyzing the vibrational energy originating from the excess contacts per particle over 
the minimum number needed for mechanical stability: Z — Zc- The extra contacts 
can be divided into Zi strong contacts and the remaining Z — Z\ weak ones. Strong 
contacts shift the energy cost of a mode that is initially zero in the isostatic limit 
to a nonzero value according to the variational argument outlined above, while weak 
contacts shift the energy simply by increasing the restoring force for displacements. 
The energy is then minimized to obtain the fraction of strong contacts. Even though 
these glasses have a high coordination number, most of the additional contacts can be 
considered to be weak. 

These results are shown in Fig. lO.lSl for two models, uj* is the theoretical prediction 
for the onset frequency of excess modes based on the variational argument. The onset 
frequency determined from simulation is given by a;^. The agreement between ui* and 
bj^ is very good. On the right-hand axis bZ = (Z\ — Zc)/Zc is shown. Note that even 
though Z itself can be arbitrarily large, 5Z remains small and does not get larger that 
0.6. This is why these systems with high coordination can still be understood in terms 
of the physics of point J and isostaticity. Even though it may be impossible to reach 
the jamming transition itself for example because there are long range interactions 
or attraction (as in the Lennard- Jones system) the effects of the long-ranged part of 
the potential can be treated as a correction to the predominant effects of jamming. 

Thermal conductivity: The thermal conductivity, k, of a system near the jamming 
threshold has been calculated and captures well the thermal conductivity of silica glass 
above its plateau( pCu et al, 2009| [Vitelli et al, 2010] ) . Because the phonons are scat- 
tered so strongly by the disorder, phonon-phonon scattering is relatively unimportant 
so the thermal transport can be calculated within the harmonic approximation with- 
out recourse to the anharmonic properties of the modes. This is in strong contrast 
to the situation for crystals, where anharmonicity must be included to obtain a ther- 
mal conductivity that is not infinite. At the jamming threshold, the long-wavelength 
phonons that give rise to a diverging thermal conductivity in the harmonic limit are 
completely suppressed by the excess modes, since the density of vibrational states re- 
mains fiat down to zero frequency. As a result, the thermal conductivity is finite even 
within the harmonic approximation at this special point. The thermal conductivity 
can be expressed in terms of the diffusivity, d{uj) (defined in section FO. 2. 6|) . and heat 
capacity, C(a;), of each mode: 

/•OC 

k{T) = I D{u) C{u)d{uj) du , (0.10) 
Jo 

In a weakly scattering system such as a crystal, d{uj) = c£{uj) /3 where c is the sound 
speed and l{uj) is the phonon mean-free path. As discussed in section I0.2.6[ Fig. 10.81 
shows that at the jamming transition d{uj) has a very simple shape: it is a plateau that 
extends from our lowest frequency up to the onset of the high-frequency localized states 
dXu et al, 2009l|Viteni et al, 2010D , where d{iv) falls rapidly to zero. The magnitude 
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Fig. 0.15 Characteristic frequencies and fractional deviation from isostaticity for two tiiree- 
-dimensional systems based on Lennard- Jones potentials. Left panel shows data for repulsive 
Lennard- Jones bidisperse mixtures where only the repulsive part of the potential is included 
versus A^, distance from the unjamming transition. Right panel shows data versus density, p, 
for a full Lennard- Jones potential with cutoff at r = 2.5aij where aij is the average diameter 
of particles i and j. On the left axis, uj^ (open circles) is the onset frequency of the anoma- 
lous modes calculated from numerical simulations and the corresponding predictions for uj* 
(crosses) based on a variational calculation. The agreement is very good. Lines are to guide 
the eye. On the right axis is the fractional deviation from isostaticity: SZ = {Z\ — Zc)/Zc 
(solid triangles). Figure taken from ( |Xu et al, 2007[ ) 



^ where 



of the diffusivity in the plateau region, do, is small and scales as: do ^ c 

keff = ^ X-^a''^^ is the effeetive spring constant of the system that scales as the bulk 
modulus, B. Effective medium theory suggests that these results hold more generally 
for weakly-coordinated disordered elastic networks ( jWyart, 2010[ ). 

Upon compression, the plateau region in ditj) no longer extends down to zero fre- 
quency but only to a frequency that is proportional to uj* . Thus all the modes above uj* 
(excluding the truly localized modes at high frequency) have a nearly constant diffusiv- 
ity. As we emphasized in the section on low-temperature excitations, uj* obtained from 
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the density of states corresponds to the frequency of the Boson peak. Here we find that 
this same frequency also is proportional to the onset of a region of flat difFusivity for 
thermal transport( Xu et al., 2009[ [Vitelli et al., 2010" |. It is just such a region of con- 
stant difFusivity that was postulated by Kittel(Kittel, 1949 1 to explain why the thermal 
conductivity of glasses has a weak, nearly linear, temperature dependence above the 
thermal-conductivity plateau. This is consistent with other evidence in glasses that 
the end of the plateau in k(T) corresponds to the Boson peak( Ze"ller and Pohl, 197T| ). 

Failure under applied stress: As summarized in section l0.2.6[ the vibrations of soft- 
sphere packings are quasi- localized near the onset of anomalous modes, uj*: they have 
large displacements of relatively few particles while the rest of the mode has a very 
small amplitude. As the sample is compressed or sheared, it will eventually go un- 
stable and rearrange into a new configuration. At zero temperature, this instability 
is governed by some quasi-localized mode in the system that goes "soft so that its 
frequency is pushed down until it reaches zero (Maloney and Lemaitre, 2006). At this 
point, the sample moves into a new configuration. One cannot predict too much from 
the purely harmonic properties of a mode about what it will do when it is pushed to 
such a large amplitude that is goes unstable. However, at the very least one can say 
that at the point of instigation, the instability will initially have a very low partic- 
ipation ratio. This is consistent with the observation that when a glass fails due to 
pressure or shear stress, the failure tends to be highly localized in shear-transformation 
zones ( |Falk and Langer, 1998| [Langer, 2008[ ). So far, however, it has not been possible 
to identify shear-transformation zones unless they flip (unless a rearrangement occurs) . 
The existence of many low-frequency quasilocalized modes with low energy barriers 
suggests that it may be possible to identify multiple shear-transformation zones from 
structural properties of the quasilocalized modes, and to predict which ones will flip 
a priori, even in large systems where there are many plane-wave-like modes at low 
frequencies. 

Structure: At point J, the pair distribution function has many singular properties. 
In particular, we have shown that the flrst peak in g(r) is a (5-function at the nearest 
neighbor distance, a. The area under this (5-function is just the isostatic coordination 
number: Z = 6 in three dimensions. On the high side of this (5-function, g(r) has a 
power-law decay: g{r) oc (r — ct)^°-^. Also, at the jamming transition, the second peak 
of g{r) is split into what appear to be two singular sub-peaks one at r = \/3a and 
the other at r = 2a. There is a divergent slope just below each subpeak and a step- 
function cutoff on its high side. The existence of a split second peak in g{r) is a feature 
that is seen in many glasses, such as metallic and colloidal glasses. One can thus trace 
the origin of these subpeaks to the geometry at the jamming threshold. One can also 
ask why, as a supercooled liquid is cooled into the glassy state, there is no signature 
in the static structure factor, S{q), which is how the structure of molecular glasses is 
determined from scattering experiments. At the jamming threshold, there are indeed 
clear signatures in the structure of a sample: these are the divergences just mentioned 
in g{r). However, because S(q) is the Fourier transform of g{r), the divergences simply 
appear as slightly more pronounced oscillations in S{q) out to very large wavevector. 
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0.6 Connection with Super-cooled liquids 

The understanding of the glass transition remains one of the enduring mysteries 
of condensed matter physics. There is no doubt that the time scale for relaxation 
increases dramatically (faster than an Arrhenius law for what are called "fragile" 
glasses) as the temperature is lowered from the melt towards the glass transition tem- 
perature, Tg, but it is much less clear whether there are any static growing length 



scales that can be identified that accompany the slowing down pkiiger et ai, 1996 
[Debenedetti and Stillinger, 2001[ ). However, as is well-documented in this book, there 
is a growing length scale associated with heterogeneities in the dynamics, which become 
increasingly collective as the viscosity increases. The origin of the super- Arrhenius in- 
crease of relaxation time with decreasing temperature (i.e. fragility) is controversial. 
Is it due to free- volume effects, thermal activation over energy barriers that grow with 
decreasing temperature, or some other type of cooperative motion? 



0.6.1 Glass transition and soft modes in hard sphere liquids 

Goldstein ( [Goldstein, 1969[ ) has proposed that the slowing down of the dynamics is 
associated with the emergence of meta-stable states as a liquid is cooled. Such a rough- 
ening of the energy landscape is predicted by mean-field models of liquids and spins 
([Mezard and Parisi, 1996 I Kurchan and Laloux, 1996 Castellani and Cavagna, 2005 



[Biroli and Bouchaud, 2009D at the so-called mode-coupling temperature Tmct- How- 
ever these theoretical approaches lead to a (non-observed) power-law divergence of the 
viscosity at Tmct, and their interpretation in real space is unclear and actively stud- 
ied ( [Franz and Montanari, 2007| [Biroli et ai, 2006| ) . Empirically, a direct validation of 



Goldstein scenario in liquids is difficult. Several numerical analysis! Angelani et ai, 2000 



Broderix et ai, 2000 Grigera et ai, 2002 1 have claimed to observe a transition in the 



energy landscape at Tmct, but in our opinion the interpretation of these results is 
disputable 0. Another intriguing aspect of Goldstein proposal is its application to hard 
spheres, where free volume rather than energy matters. 

The connection between microscopic structure and vibrational spectrum estab- 
lished above directly supports that the emergence of meta-stable states slows down the 
dynamics, and yields a geometrical interpretation of this phenomenon ( [Brito and Wyart^ 
In order to see that, the theoretical description of anomalous modes and their as- 
sociated length scales must be extended to finite temperature. An analogy between 
hard-sphere liquids and athermal elastic systems can be built ( Brito and Wyart, 2006[ 



Brito and Wyart, 2009) in two steps: (i) on short time scales particles rattle rapidly 



2009] ). 



^ In these numerics, equilibrium configurations close to the configurations visited by the dynamics 
are analyzed. Their number of saddles (or unstable modes) is computed in function of the temperature, 
and it is argued that this quantity can be fitted by a curve that vanishes at a temperature Tmct, 
defined here by fitting the dynamics as r ~ {T — Tmct)'"" , where t is the a-relaxation time scale and 
a some exponent. From this observation it is stated that meta-stable states appear at Tj^jct- However, 
this observation is precisely the behavior expected if meta-stable states had already appeared at larger 
temperature, implying that the dynamics is activated in the neighborhood of T^CT- Indeed assuming 
activated dynamics, the number of saddles ris reflects the probability that a region is in the process 
of rearranging (crossing a barrier), which occurs with a probability of order 1/t. Thus one expects 
ris ~ 1/r, which can be fitted with a vanishing quantity at Tj^jCT simply due to the definition of 
that temperature. 
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Fig. 0.16 Left: Illustration of the fact that the problem of hard elastic discs at finite temper- 
ature and density in the supercooled liquid or glass phase (in this case 4>c — 4> ~ 2-10^'' 0c), can, 
in between rearrangement events, be mapped onto the problem of a static packing with a log- 
arithmic interaction potential ( [Brito and Wyart, 2006| |Brito and Wyart, 2007[ ) . Essentially, 
particles experience a reduced free volume due to collisions with close neighbors, leading to an 
effective repulsive interaction. Thus particles exert a force on each other. The plot shows the 
average contact forces for 256 particles averaged over 10^ collisions. The points represent the 
centers of the particles, and the thickness of the line is proportional to the force strength. Since 
in between rearrangement events the structure is stable, forces are balanced on each particle 
( |Brito and Wyart, 2006 1 [Brito and Wyart, 2007p . Right: IS.Z inferred from the coordination 
of the force network vs. average contact force (/) , a measure of pressure. Measurements are 
made in the supercooled liquid (diamonds) and in the glass (circles). The observed coordi- 
nation is well-captured by its minimal value AZ ~ (f)^^^^ allowing mechanical stability, 
supporting that the glass is marginally stable ( [Brito and Wyart, 2009 [ I. 

but no large rearrangements occur. As illustrated in Fig. 10.161 a contact force network 
can be defined, following earlier work on granular matter ( [Ferguson et ai, 2004| ), which 
represents the average force exchanged between particles. From such networks a co- 
ordination number is computed by counting the pairs of particles with a non-zero 
contact force, (ii) All the possible hard-sphere configurations associated with a given 
network can be summed up, leading to a computation of the free energy that can be 
expressed in terms of the mean positions of the particles. The effective interaction is 
found to be logarithmic. Expanding the free energy defines a dynamical matrix and 
vibrational modes that characterize the linear response of the mean particles positions 
to any imposed external forces within a meta-stable state. 

The results on the vibrational spectra of athermal elastic networks apply to this 
case as well, and the stability criterion of Eq. (|0.8p must be satisfied in any metastable 
configuration. For hard spheres, it can be written AZ > p~^^^ ^ h^^^ as the pressure 
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Fig. 0.17 Phase diagram for tlie stability of hard-sphere configurations, in the coordination 
5z vs average gap h plane. The marginal stability line delimits stable and unstable con- 
figurations. The dashed line correspond to equilibrium configurations for different (j). As </> 
increases, h decreases and the two lines eventually meet. This occurs at the onset packing 
fraction 4'onset, where the dynamics become activated and thus intermittent. At larger 0, 
viscosity increases sharply as configurations visited become more stable. For a finite quench 
rate the system eventually falls out of equilibrium at the glass transition packing fraction 
00 ■ More stable and more coordinated regions cannot be reached dynamically, and as is 
increased further, the system lives close to the marginal stability region, as indicated in the 
dotted line. The location of the out-of-equilibrium trajectory depends on the quench rate. In 
the limit of very rapid quench, the out-of-equilibrium line approaches the marginal stability 
line dBrito and Wyart, 2009"| . 

p satisfies p ^ 1/h ^ where h is the typical average gap among particles 

in contact. As illustrated in Fig. 10.161 the bound is approximatively saturated and 
the glass appears to have just enough contacts to maintain its stability. Marginal 
stability is a fundamental feature of the glass, generically foreign to crystals. It also 



occurs in some mean field spin models (Kurchan and Laloux, 1996). In hard spheres it 



implies the presence of anomalous modes near zero frequency and a curious behavior 
of the short-time dynamics, in particular a dependence on the mean-square particle 
displacements scaling as {SR"^) ~ Ac/)^/^ rather than the AcjP dependence expected in 



a crystal (Brito and Wyart, 2009). Those predictions have been confirmed numerically 
dBrito and Wyart, 2007||Brito and Wyart, 2009l|Mari et al, 2009D. 



Why is a hard sphere glass marginal? It has been argued ( Brito and Wyart, 2009[ ) 
that this must be so if the viscosity increases very rapidly when meta-stable states ap- 
pear in the free-energy landscape. In the (AZ, h) plane sketched in Fig. 10.171 Eq. (|0.8p 
draws a line separating stable and unstable configurations. At equilibrium the mean 
value of AZ and h arc well-defined functions of 0, and the corresponding curve is 
represented by the dashed line in Fig. 10.171 At low 0, gaps between particles are large 
and the configurations visited are unstable. As increases, the gaps decrease and con- 
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Fig. 0.18 Above: two examples of sudden rearrangements in the glass phase for different 
average contact force (/) for A'^ = 1024 particles. Below: projection of these rearrangements 
on the 1% of the normal modes that contribute most to them, which systematically lie at the 
lowest frequencies ( |Brito and Wyart, 2009p . This projection captures most of the displace- 
ment, indicating that a limited number of (collective) degrees of freedom are active during 
rearrangements. 

figurations eventually beeome stable. This occurs at some (ponset when the equilibrium 
line crosses the marginal-stability line. At larger 0, the dynamics becomes activated 
and intermittent, and the viscosity increases rapidly according to our hypothesis, so 
that on accessible time scales equilibrium cannot be reached deep in the regions where 
metastable states are present. As a consequence, the system falls out of equilibrium 
at some 0o larger than (ponset- Configurations visited must therefore lie close to the 
marginal stability line, as represented by the dotted line in Fig. 10.171 since more stable, 
better-coordinated configurations cannot be reached dynamically. 

Further support for the view that the viscosity is governed by the presence of 
an elastic instability near (j)onset comes from the observation that the sudden rear- 
rangements that relax the structure occur along its softest modes (both in the liq- 
uid and in the glass phase) ( |Brito and Wyart, 2007[ ), indicating that the lowest ac- 
tivation barriers lie in these directions, as illustrated in Fig. 10.181 Experiments in 



granular materials (Brito et al., 20101 and numerical observations of normal modes of 



the energy of inherent structures ( Schober et al., 1993 Widmer-Cooper et al., 2009a 



[Widmer-Cooper et al., 2009fe[ ) reach similar conclusions. The normal modes of the free 
energy considered here are expected to correlate with the dynamics much better than 
those of the energy, and must be considered in particular when the interaction po- 
tential is strongly non-linear. The collective rearrangements involving a few tens of 
particles commonly seen in super-cooled liquids correspond mostly to a few modes 
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Fig. 0.19 Relaxation time, r, versus temperature, T, (a) and versus inverse pressure, 1/p, (b) 
for a system of spiieres interacting via repulsive harmonic potentials. Panel (c) shows that all 
the relaxation time data for different potentials (harmonic, Elertzian and hard sphere) can be 
collapsed onto a single master curve. To get this data collapse, the relaxation time is non-di- 
mensionalized by pressure and plotted versus T/pa^. Figure taken from ( |Xu et ai, 2009[ ). 



dBrito and Wyart, 20071 ), an effect even stronger in the glass phase. This analysis sup- 
ports that the anomalous modes are the elementary objects to consider to describe 
activation, and that the collective aspect of rapid rearrangements is already present 
at the linear level in the structure of the modes. 

This approach permits a quantitative study of activation, in particular to com- 
pute the fraction of the spectrum contributing to rearrangements. This quantity de- 



creases rapidly with compression near the glass transition (Brito and Wyart, 2007), 
implying that fewer and fewer degrees of freedom are relevant for the dynamics. This 
observed reduction of active degrees of freedom presumably mirrors the slowing down 
and increasingly collective aspect of the rapid rearrangements that relax the structure 
( [Candelier et ai, 2010[ |Candelier et ai, 2009D , but this has yet to be clarified. 

Overall, the present approach furnishes a simple spatial picture for the slowing- 
down of the dynamics at intermediate viscosities: meta-stable states appear when the 
contact force network become sufficiently coordinated to resist the destabilizing effect 
of pressure. The localization for long times of an individual particle is thus not due to 
some specific properties of the cage formed by its neighbors, but has to do with the 
structure of the packing on a length I* that can be large. Incidentally the concept of 
caging, commonly used to picture dynamical arrest, is misleading. 



0.6.2 Models of soft finite range repulsions at finite temperature 

By analyzing several models near the jamming threshold, some universal aspects of the 
slow-relaxation dynamics that appear in one asymptotic limit become apparent ( |Xu et ai, 2009 ). 
In particular, when the pressure, p, is small all the relaxation time data can be suitably 
plotted on a single master curve as a function of a single variable. 

The models considered here all have finite-range repulsions. For this study, thermal 
energy was included in the simulations, and the relaxation time r of the system was 
measured by studying when the intermediate scattering function, f{q,t) drops to 
of its initial value. 

Normally in simulations of this type, one measures time in units of a^/rn/e where 
a is the particle diameter, m is the mass of the particle, and e is the energy scale 
of the potential. However, in order to produce a collapse onto a master curve, r is 
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non-dimensionalized in another way by using Ty^pa/m. Plotting this quantity versus 
T/{p<T^) produces an excellent collapse of the data for several different potentials of 
interaction at all values of T/ (pcr^) as long as the non-dimensionalized pressure, pa^ /e 
is sufficiently small( Xia et ai, 2009[ ). This is shown in Fig. 10.191 This limit is satisfied 
at low pressures in soft-sphere systems and at any pressure for a hard-sphere system 
where e — > oo. Fig. 10.191 shows that low-pressure soft-sphere data do indeed collapse 
with the hard-sphere result. In addition, Ref. ( |Xu et ai, 2009[ ) found that limiting 
hard-sphere behavior remains a good approximation up to higher values of ptr^/e as 
T/{p(T^) increases. 

This data collapse has two important implications. First, it shows that at suf- 
ficiently low pressures the relaxation at the glass transition is a function of only a 
single variable. The molecular glass transition is, in this limit, equivalent to the hard- 
sphere colloidal glass transition. Second, Ref. (Xu et ai, 20091 shows that there are 
systematic deviations from this universal behavior at larger values of the second con- 
trol parameter: pa^ /e. This demonstrates that there are at least two distinct physical 
processes that enter to produce glassy dynamics. What is important is that they can 
be cleanly separated, so that for sufficiently small pa^ /e, the rescaled relaxation time 
is only a function of a single variable, T/{pa^). Note that the existence of such a 
scaling function still does not tell us whether the rescaled relaxation time diverges at 
T I (pa^) > 0, implying that a thermodynamic glass transition exists, or at T/ (pa^) = 0, 
implying that the glass transition is an isostatic jamming transition. 

In order to understand the consequences of this result for molecular liquids, we 
must first understand the corrections to the leading hard-sphere behavior when we 
slowly increase pa^/e. These studies, have been done for repulsive spheres near the 
jamming threshold. However, molecular liquids typically have attractions and are at 
much higher densities than those considered here. To make the correspondence more 
general, we would have to push these systems farther from the jamming threshold. 
Nevertheless, this result shows how certain aspects of jamming, inherent in the hard- 
sphere liquid, are important for glass-forming liquids generally. As we have shown, at 
least one other distinct contribution to the relaxation must also be considered as the 
pressure is increased. 

As one increases pa'^/e, there are corrections to the leading hard-sphere behavior. It 
has been argued that these corrections can be collapsed onto a single curve by rescaling 
the pressure by some factor that varies with packing fraction ( [Berthier and Witten, 2009 6[ ) . 



0.7 Outlook to the future - A unifying concept 

The vibrational spectrum of crystals consists of plane waves and leads to a descrip- 
tion in terms of linear elasticity. In those systems, one can compute the non-linearities 
of the excitations and how they scatter from defects. Such a systematic description 
is lacking in amorphous solids. Various properties, such as thermal transport, force 
propagation, the vibrational spectrum and resistance to flow present phenomena that 
are not satisfactorily understood. An inherent difficulty in understanding these phe- 
nomena is that they are often governed by processes that occur on small length scales, 
of the order of a few particle sizes. 
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The discovery that the jamming threshold of athermal ideahzed spheres exhibits 
many properties of a critical point ( [O'Hern et ai, 2003] ) changes this state of affairs, as 
it enables one to separate the length scale on which disorder matters from the molec- 
ular scale. This has led to the realization that the spectrum of amorphous solids is 
characterized by a crossover, distinct from Anderson localization, above which modes 
are extended but not planc-wave-like, and that the characteristic length at which 
this crossover occurs is decoupled from the molecular length. The corresponding ex- 
citations, the anomalous modes, allow one to understand several seemingly disparate 
anomalies in glasses and granular materials in a coherent fashion. The picture ad- 
vanced in this review suggests that (1) the length scale I* corresponding to the 
crossover characterizes force propagation in granular packings ( [Ellenbroek et ai, 2006[ 
|Ellenbroek et ai, 2009 a[ ), (2) the frequency scale of the crossover corresponds to that 
of the boson peak in glasses ( [Silbert et ai, 200"5 Wyart et ai, 2005a Xu et ai, 2007 



|Wyart, 2005[ ), (3) the low-energy diffusivity of the modes above the crossover is re- 
sponsible for the mild linear increase of the thermal conductivity with temperature 
above the plateau in glasses ( |Xu et ai, 2009] |ViteUi et ai, 2010| [Wyart, 2010D , (4) at 
least for colloidal glasses, the microscopic structure is marginally stable toward anoma- 
lous modes ( [Brito and Wyart, 2006|[Brito and Wyart, 2009[ ), (5) the modes appear to 
govern structural rearrangements in liquids at least in the viscosity range that can be 
probed numerically ( [Brito and Wyart, 2007 Brito and Wyart, 2009 Widmer-Cooper et ai, 2009a land 
(6) the structure of the modes corresponds to that of structural arrangements in amor- 
phous solids under mechanical load ( [Xu et ai, 2010| [Manning and Liu, 2009[ ). At the 
linear level those excitations and their consequences are now becoming rather well- 
understood theoretically for sphere packings, with the notable exception of (i) the 



observed quasi- localization of the lowest- frequency anomalous modes (Xu et ai, 2010 



[Zeravcic et ai, 2d09[ ) and (ii) the apparently fractal statistics of force propagation 
below the cut-off length £* above which continuum elasticity applies. 

The jamming community is now beginning to study the non-linearities coupling 
these excitations. Much would be gained from such knowledge. On the one hand, it is 
possible that the lowest-frequency anomalous modes are related to two-level systems, 
which have been proposed on phenomenological grounds to explain transport and 
specific heat anomalies at sub-Kelvin temperatures. Despite decades of research, the 
spatial nature of two-level systems has remained elusive. If there is indeed a connec- 
tion between the highly anharmonic, low-frequency, quasilocalizcd modes and two-level 
systems, then the energy barrier separating the two levels should vanish as the iso- 
static jamming transition is approached. On the other hand, since the low-frequency 
modes are involved in irreversible rearrangements, understanding their non-linearity 
is presumably essential for a microscopic description of the flow of supercooled liquids, 
granular matter and foams. 

An exciting recent development is the advent of experiments capable of studying the 
vibrational properties of systems near the jamming transition ( [Abate and Durian, 2006 



[Ghosh et al, 2010[ [Chen et al, 2010[ [Gardel, 2009t . In laboratory systems, particle 
motion is typically damped and the inter-particle potentials are not necessarily known; 
this has made it difficult to make direct comparisons between experimental results and 
numerical and theoretical predictions until now. These complications have recently 
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been overcome ( |Ghosh et ai, 2010] |Chen et ai, 2010| , so experiments should soon be 
able to test the extent to which predictions for ideal spheres apply to real systems. 
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